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Abstract. We develop a Non-Crossing Approximation (NCA) for the effective cluster problem of the re- 
cently developed Dynamical Cluster Approximation (DCA). The DCA technique includes short-ranged 
correlations by mapping the lattice problem onto a self-consistently embedded periodic cluster of size Nc- 
It is a fully causal and systematic approximation to the full lattice problem, with corrections C'(l/A'^c) in two 
dimensions. The NCA we develop is a systematic approximation with corrections 0{1/N^). The method 
will be discussed in detail and results for the one-particle properties of the Hubbard model are shown. Near 
half filling, the spectra display pronounced features including a pseudogap and non-Fermi-liquid behavior 
due to short-ranged antiferromagnetic correlations. 

PACS. 7 l.lO.Fd - 7 1.27.-Fa - 7 5.20.Hr - 7 5.30.Kz - 7 5.30.Mb 



1 Introduction 

One of the most challenging tasks in theoretical condensed 
matter physics is the description of strongly correlated 
electron systems. The Coulomb interaction between elec- 
trons plays a dominant role in these systems and usu- 
ally strongly influences the electronic properties and the 
physics of the ground-state phase. The discovery of heavy 
Fermion systems in the late seventies [El and high-Tc super- 
conductors ten years later [|| has stimulated strong ex- 
perimental and theoretical interest in this field. However, 
despite a multitude of attempts to describe such strongly 
correlated electron systems theoretically, a complete un- 
derstanding of the observed rich physics has not yet been 
accomplished. Even the simplest model for strongly cor- 
related electron systems, the Hubbard model (HM), must 
be considered unsolved in more than one dimension [BJ af- 
ter almost forty years of intensive study. Exact diagonal- 
ization or quantum Monte-Carlo studies for two or three 
dimensions are restricted to small lattice sizes and predic- 
tions for the thermodynamic limit are maybe problematic 

i- 

However, in the limit of infinite dimensions D = oo, 
correlated lattice models undergo a significant simplifica- 
tion. Their dynamics become purely local and therefore 
the lattice problem can be mapped onto a generalized sin- 
gle impurity Anderson model coupled to a host, which 
has to be determined self-consistently P,0M@1. The Dy- 
namical Mean Field Approximation used in the context 
of real materials thus assumes that only local dynamics 



are present. Despite the neglect of nonlocal correlations, 
this method has been shown to capture several important 
features of e.g. the Hubbard model M,0,M]. Nevertheless 
it has some significant shortcomings due to the mapping 
onto a purely local model. For instance, it does not in- 
clude the effect of nonlocal correlations like antiferromag- 
netic spin fiuctuations on the one-particle properties, and 
is not capable of describing nonlocal order parameters. 
However, both effects are believed to be especially im- 
portant for a description of the High-Tc materials. Here 
the one-particle spectra have shadow bands due to short- 
ranged antiferromagnetic fluctuations, preformed pseudo- 
gaps due to superconducting or spin fluctuations [0,Q3, 
0,n^na|, and the superconducting order parameter is of 
nonlocal (d) character. 

In order to include these types of nonlocal dynamics 
into the theory there have been several efforts to add so- 
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called -^ corrections to the DMFA |]J,|l7|,|l8|. However, 
these methods either experience causality problems [ |l6| , 
|l7| (because of the necessary inclusion of nonlocal Green 
functions in self-energy diagrams which do not have a neg- 
ative semidefinite imaginary part), or are restricted to the 
calculation of a few moments of the spectral function [|9[ . 
These shortcomings do not apply to the Dynamical 
Cluster Approximation (DCA). This method systemati- 
cally includes nonlocal short-ranged correlations while pre- 
serving causality |g^,|l|. The DCA is a scheme which 
maps the lattice problem onto a self-consistently embed- 
ded effective finite-size cluster model. Due to the finite 
size of the cluster, nonlocal corrections to the local dy- 
namics can be systematically included as the cluster size 
increases. The basic idea of the DCA is to take into ac- 
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count nonlocal physics by calculating the self energy at 
selected points K in the Brillouin zone and consider the 
self energy at these points to represent the self-energy in 
the surrounding of these points K -I- k': S(K) ~ ^(K-|-k'). 
The theory then maps the lattice problem onto an effec- 
tive finite-size system with periodic boundary conditions 
coupled to an external bath and the resulting system is 
solved self-consistently. The DCA has two well defined 
limits: ft recovers the DMFA as the cluster size goes to 
1 and becomes the exact solution for the model under 
consideration as the cluster size goes to infinity. 

The use of the DCA as an approximation can be jus- 
tified as long as the momentum dependence of the self- 
energy of the real system is weak. This is obviously re- 
alized in high spatial dimensions where a coarse grid of 
K-points should capture all the basically short-ranged dy- 
namics. In two- or three-dimensional systems the approx- 
imation is more crude, but can be motivated by the ob- 
servation that the dominant structures in the one-particle 
dynamics are generated by local renormalizations, while 
nonlocal effects only lead to minor renormalizations of 
these structures. Note that this assumption automatically 
inhibits studies very close to phase transitions since there 
strong, long-ranged fluctuations must be expected. How- 
ever, sufficiently away from phase boundaries correlated 
systems indeed show only mild momentum dependence of 
the one-particle self energy as compared to its frequency 
dependence. 

So far only Quantum Monte Carlo simulations and ex- 
act enumeration have been used to solve this problem of 
a cluster in an external bath [|lH|§. In the DMFA the 
NCA has successfully been applied to the effective single 
impurity problem ^,^,24,2^,^. In this paper we intro- 
duce an extended version of the NCA to solve the effective 
periodic cluster model of the DCA. 

The paper is organized as follows. First a short review 
of the DMFA is given, which is reproduced by the DCA for 
a single site cluster. Then we provide a microscopic defini- 
tion of the DCA in terms of its Laue function, and rederive 
the DCA algorithm using Baym's <? functional formalism. 
We then define the effective cluster model onto which the 
lattice system is mapped by the DCA. An extended ver- 
sion of the NCA applicable to an impurity cluster of arbi- 
trary size is discussed in detail and finally results for the 
one-particle properties of the Hubbard model are shown 
and compared to corresponding results of the DMFA. 



2 Dynamical Mean Field Approximation 

In this paper we consider the single-band Hubbard model 
described by the Hamiltonian 



in infinite dimensions, while retaining their full local dy- 
namics. Mctzncr and Vollhardt|g^ showed that the nec- 
essary rescaling of the kinetic energy as 1/vD leads to 
a collapse of all nonlocal diagrams in a skeleton expan- 
sion for the self-energy. Consequently the corresponding 
Baym-Kadanoff <P functional can be expressed in terms of 
local quantities only. 

Miiller-Hartmann p8[ was able to deduce the same 
result by inspecting the momentum dependence of vertices 
in diagrammatic approaches as D — *■ oo. For Hubbard-like 
models, the momentum dependence of each vertex in a 
diagrammatic expansion of the functional <? is completely 
characterized by the Laue function 



Z\(ki,k2,k3,k4)=^e^('^^-''=+'^-^ 



(2) 



where ki and ka (k2 and k4) are the momenta enter- 
ing (leaving) the vertex. In a conventional diagrammatic 
approach Z\(ki,k2,k3, k4) = A^5ki+k3,k2-i-k4, which ex- 
presses momentum conservation on the vertex. However 
as £) — > oo Miiller-Hartmann showed that the Laue func- 
tion reduces to 



Z\,3^oo(kl, k2, k3, k4) = 1 + 0{1/D) 



(3) 



The DMFA assumes the same Laue function (H) even 
in the context of finite dimensions. Therefore both the 
infinite-dimensional theory and the DMFA neglect mo- 
mentum conservation at the internal vertices of irreducible 
diagrams and the momenta in the corresponding 'Pdmfa 
functional may be freely summed over the whole Brillouin 
zone. This leads to a collapse of the momentum depen- 
dent contributions to the functional 'Pdmfa and only lo- 
cal terms remain. This is illustrated in Fig. || for a second 
order diagram. 



k4 
k. 
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A = A 



(i=oc 



ki 






Fig. 1. The Laue function ra) leads to a collapse of the momen- 
tum dependence of the DMFA Baym-Kadanoff '^dmfa func- 
tional, illustrated for a second order contribution. 



H^Y.^^AaCJa+UY.n,^ 



"a 



(1) 



where c\^ (cj^^) creates (destroys) an electron at site i with 
spin sigma and riifj are the corresponding number oper- 
ators. Lattice models of this kind simplify significantly 



The self-energy (a functional derivative of the func- 
tional ^DMFA with respect to a Green function leg) also 
becomes a functional of local propagators only and there- 
fore becomes a constant in momentum space. Consequently 
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the lattice problem can be mapped onto an effective im- 
purity problem. 

The DMFA has proven to capture the key features of 
strongly correlated electron systems and to provide insight 
in the complicated dynamics mediated by correlations. De- 
spite its great success in the description of correlated elec- 
tron systems the DMFA has some significant shortcomings 
due to the neglect of non-local dynamics. 



3 Dynamical Cluster Approximation 

The DCA extends the DMFA through the inclusion of 
short-ranged dynamical correlations. The DCA was intro- 
duced and discussed in detail in [g0||2^. In this paper we 
will rederive the DCA algorithm with an argument which 
is complimentary to that used by Miiller-Hartmann |g8[ 
to describe the DMFA. 
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of the Laue function the momenta of each internal leg in 
the corresponding functional 'Pdca niay be freely summed 
over the cell and each leg is replaced by the coarse grained 
average 



k 



k) 



(5) 



This is schematically illustrated in Fig. ||. 
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Fig. 3. The DCA choice of the Laue function (Ul) leads to 
the replacement of the lattice propagators G(ki), Gfk2), ... 
by coarse grained propagators G(K), G(K'), ... (Eq. |5|) in the 
internal legs of <Pdca, illustrated for a second order diagram. 



The coarse grained Green function G(K) and correspond- 
ing self-energy S(K.) 



N, 



2n/L 



^(K) = ^E^(K+k) 



(6) 



Fig. 2. The cluster momenta K (filled circles) and coarse 
graining cells (different fill patterns) for a A'^c = 4 cluster in 
the Brillouin zone of a two dimensional lattice. The cells adja- 
cent to the Brillouin zone boundary extend periodically to the 
opposite site. 



The basic idea of the DCA is to partially restore the 
momentum conservation relinquished by the DMFA. To 
this end the Brillouin-zone is divided into Nc — L^ cells 
of size {2'k/L)^ (see Fig. g). Each cell is represented by 
a cluster momentum K in the center of the cell. We re- 
quire that momentum conservation is (partially) observed 
for momentum transfers between cells, i.e. for momentum 
transfers larger than Ak ~ 27r/L, but neglected for mo- 
mentum transfers within a cell, i.e less than Ak. This re- 
juirement can be established by using the Laue function 



quiri 

m 



^DCA 



NrS- 



M(ki)+M(k3),M(k2)+M(k4) 



(4) 



where ]V[(k) is a function which maps k onto the cluster 
momentum K of the cell containing k. With this choice 



are functions of the cluster momenta K only. The k sum- 
mation in (P) and (m runs over the N/Nc momenta of the 
cell about the cluster momentum K and G (E) is the full 
lattice propagator (self-energy). 

With this definition of Adca, the DCA estimate of 
the Baym-Kadanoff functional 'Pdca becomes pll 



1 



^I3Ca(G) -^-tr ^«G 



(7) 



where Z'^'^ is the set of irreducible, coarse grained self- 
energy diagrams of l^^ order in the interaction U and 
the trace indicates summation over frequency, cluster mo- 
menta and spin. The DCA result for the free energy of the 
lattice is 



nocA = ~kBT{$DCA - tri;G- trln[-G]) 



(8) 



^DCA is stationary with respect to the lattice Green func- 
tion G 

^ = fc.Thi: + r] = o , (9) 

if S is taken as an approximation for the self-energy S of 
the full lattice Green function G (the left hand side of (||) 
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follows from SCk/SG^ = <5K,M(k) and S<Pdca/^G = S). 
We have shown previously that S — S + 0{l/Nc) in two 
dimensions and includes dynamical intersite correlations 
of range 7r/Zifc = L/2 |l|. 

The coarse grained Green function ^ then takes the 
form 



Nr. 



'5(K,.) = ^E 



^K+k 



/i - S{K, z) 



(10) 



where the self-energy at momentum k = K + k, Z'(k) 
is replaced by its coarse grained average S{K.) and z = 
uj + iS. Note that the choice of the coarse grained Green 
function (|o|) has two well defined limits with respect to 
the cluster size Nc- For Nc = 1 the k summation runs 
over the entire first Brillouin-zone, G is the local Green 
function, thus the DMFA algorithm is recovered. When 
Nc = oo the k summation vanishes and the DCA becomes 
equivalent to the exact solution of the Hubbard model. 

In order to apply the NCA to solve the effective clus- 
ter problem it is convenient to write the coarse grained 
Green function ([lOl ) in a more suitable form. We use the 
independence of the self-energy i?(K) on the integration 
variable k to write G in the form 



G'(K,z) 



1 



z - CK + M - '£'(K, z) - r(K, z) 



(11) 



where Iml^(K) < (see Appendix ^). This is just the 
Green function of an effective cluster model with peri- 
odic boundary conditions coupled to a Fermionic bath de- 
scribed by the host function r{K.). Hence we can obtain 
the coarse grained self energies ^CK.) by solving a gener- 
alized cluster model. 

The DCA cluster problem may then be solved by iter- 
ation. The iteration loop starts with a guess for the initial 
cluster self-energy. By computing the coarse grained lat- 
tice Green function ( |lO| ) we get the input quantities for 
the effective cluster model. The effective dispersion in is 
given by the average 



N, 



EK = 



k 



(12) 



and the cluster electrons are coupled to the host func- 
tion -r(K). This will be described in more detail in sec- 
tion and Given the effective dispersion ck and the 
host r(K) the interacting Green function G'c(K) of the 
effective cluster model can be calculated by some suit- 
able method. The cluster self energy is then obtained via 
Ec{K,z) = z - ex + A* - r{K,z) ~ G-^{K,z) and the 
iteration closes by calculating a new G(K) with Eq. HO. 
This procedure is repeated until G'c(K) — G(K) within 
the desired accuracy. 



the Hamiltonian are given by the Green function (nTl) 
The corresponding cluster Hamiltonian 



H„ 



(13) 



is most conveniently constructed in momentum space. Its 
local part is given by 



Hioc - 22 '^Kfnafvia 



Kct 

U 



"*" w Z^ /k+qt/kt/k'-qx/k' 



I ' 



(14) 



Q 



where /^^ (/ko-) creates (destroys) an electron with mo- 
mentum K and spin a. U is the local Coulomb repulsion 
for two electrons residing on the same cluster site. Since 
this interaction is local, it is unchanged by the coarse- 
graining procedure [pl|. Note that the effective dispersion 
ek of the cluster states is given by the average bare dis- 
persion in the cell ([l^ ) . The coupling of the local cluster 
states with the host has the form 



Joined — 



^ E ^^MfL 



o-'-K+k'o- 



+ h.c.) 



K,k'. 



+ E '^kcLcka 



(15) 



kcr 



where cj^ i k' (^K+k' ) describe the effective medium in terms 
of free fermions with a dispersion e^. Note that in contrast 
to the single impurity model, the local states given by /^^ 
couple only to fermions with momenta k = K -I- k' within 
the cell about the cluster momentum K. Therefore the 
corresponding hybridization function 



r'{K,z)^ 



\Vk.i 



N ^ z- ex+k' 



(16) 



becomes K-dependent and the interacting cluster Green 
function finally reads 



Gc(K,z) = 



1 



z-eK + /U-^c(K,z)-r'(K,z) 



(17) 



Z'c(K, z) denotes the proper one-particle self energy ef- 
fects due to the local Coulomb repulsion between the /- 
electrons. 

Comparing this result with the coarse-grained Green 
function of the lattice ( pi] ) one finds identical structures 
provided that the cluster self-energy X'c(K) equals that of 
the lattice -£'(K), and that 



r'(K) = r(K). (18) 

4 Effective Cluster Model 

The latter substitution ensures that the solution of the 
To solve the cluster problem with the NCA, we must first effective cluster model is also the solution of the coarse- 
define a Hamiltonian for the cluster. The parameters of grained lattice problem discussed in the last section. 
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5 Extended Version of the Non Crossing 
Approximation 

In the following we show how to solve the effective cluster Vf^") 
model with an extended version of the NCA. The effective 
cluster model is defined by the Hamiltonian (|l^) with the 
effective medium fixed by Eq. |IJ. The NCA is a perturba- 
tional expansion around the molecular limit, i.e. it starts 
with the eigenstates of the local part Hioc of the cluster 
Hamiltonian ([f3[). The expansion is performed with re- 
spect to the coupling to the effective medium Hmed, where 
the quasi-free fermions are described by the host function 
r'(K), see Eq. |l^. The Fermionic operators defined on the 
cluster are expanded in terms of the Hubbard operators 
Xynn = \m) (n|, e.g. 



/ko- ^ 7 , ^mn -^n 



(19) 



where {|m.)} are the eigenstates of the local Hamiltonian 



Hi, 



7 ^ EmXn 



(20) 



with eigenenergies Em and F^^ — {m\f^^\n). The hy- 
bridization term (|l5|) becomes 



^med - -^ 2J Z_/^K,k'(cK-Hk'^mn-''^r, 



h.c.) 



+ E-^ 



K,k'(T m,n 
EkC^^Cktr 



(21) 



k(T 



Since the Hubbard-operators do not obey standard Ferm- 
ionic or Bosonic commutation relations, the conventional 
Feynman diagram technique cannot be used for a per- 
turbation expansion and the concept of resolvents must 
be introduced instead |29 . Their matrix-elements in the 
space of the local eigenstates have the form 



[R ^]mn{z) = (z- Em)5„ 



.(^) 



(22) 



In general non-diagonal elements Rrmi of the resolvent ex- 
ist, but if the hybridization term H„ied does not break the 
symmetry of the cluster local Hamiltonian Hioc they are 
zero. E describes self-energy effects due to the hybridiza- 
tion with the effective medium. Note that S collects the 
renormahzations of the individual local states {|m)} and 
must not be confused with the proper one-particle self en- 
ergy of the cluster, i^(K, z). 

In the NCA the self-energy matrix S is obtained by 
calculating the two diagrams illustrated in Fig. ^, which 
correspond to 



77" r rf 



Ko 



de /(£)Imr(K, e)F'^'^R{z + e){F'^'') 



^K<TNt 



-\-OQ 

I de /(-e)Imr(K, e){F^'')^ R{z - e)F^° 



r(K) 



r(K) 




Kff 



Fig. 4. NCA self-energy of the resolvent 7? due to the coupling 
of the local cluster states to the host -r(K). 



(23) 



The coupled singular integral equations (g2|) and (|2^) have 
to be solved self-consistently. 

Higher order corrections to these equations come in the 
form of vertex corrections or crossing diagrams. For the 
orbitally non-degenerate single impurity Anderson model 
it is well known that to obtain the correct value for the 
low-energy scale one has to sum all diagrams up to fourth 
order in the coupling V^lQ. For a cluster of size Nc > 
1 and finite value for U, this requirement would mean 
that one must include vertex corrections. From our former 
experience with vertex corrections for the single impurity 
case {Nc = 1) we expect a strong renormalization of low- 
energy scales. However, for high energy features in the 
spectra or high temperature properties like magnetism the 
vertex corrections yield only negligible effects. They also 
do not affect general local properties like universality or 
scaling jsH. 

In the present context, the hybridization strength is 
not an adjustable parameter, so it does not make sense 
to use it to classify the higher-order corrections. In fact, 
both the effective hybridization strength between the clus- 
ter and its host, and the degeneracy and magnitude of 
the cluster states depend upon Nc- Therefore, a far more 
important expansion parameter is the inverse cluster size 
l/Nc- Since the eigenenergies Em of the cluster scale as 
Nc, the resolvent behaves like R ~ 0{l/Nc). Taken to- 
gether the sum over the cluster momenta K and the re- 
solvent in (g^) are of order one. Thus the A^c-dependence 
of the NCA-self energy matrix X! is determined by the Nc- 
dependence of F. We show in Appendix |b| that the host 
function F is of order 0{l/Nc). Therefore for Nc ^ oo the 
NCA-self energies vanish - as expected since the coupling 
to the host vanishes - and the cluster problem is solved ex- 
actly by the eigenstates {|to)}. In the form (^) the NCA 
equations are exact up to the second order in the coupling 
(11), i.e. first order in T - 0{1/Nc). 

To estimate the role of vertex corrections we show in 
Fig. ^ one of the leading order corrections to the NCA- 
self energy (23). This crossing diagram involves two F- 
lines and three resolvents R, but only two sums over the 
cluster momenta K, therefore this diagram is of order 
0{l/N^). In fact all crossing diagrams are of this order or 
higher. Hence for Nc ^ oo the NCA algorithm becomes 
exact with corrections 0{1/Nc)- We are thus confident 
that at least the qualitative aspects of our results will be 
unaffected by higher order diagrams. Since on the other 
hand an inclusion of vertex corrections is associated with a 
tremendous numerical effort, we refrain from taking them 
into account for the time being. 
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pKa 
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Fig. 5. Leading order crossing diagram is of order 0(1/N, 
since it involves five quantities of order C'(l/A'^c) but only two 
sums over the cluster momenta K. 



With the expansion ( |19D the cluster Green function Gc — 

((/ko-j /ko-)) *^^^ ^^ written in terms of the Hubbard op- 
erators as 

Gc(K,z) = } ^ F„^^^F^,^,{{Xmn]Xm'n'))z ■ (24) 



Within the NCA, the correlation function on the right- 
hand side of (|2^) can be written as 

4-00 

{{Xmn;Xm>n'))uj = -^ / ds e~^^(p„'m (e)i?„m' (s -I" w) 
— oo 

- Pnm'ie)Rn'm{£ -Uj)) ■ (25) 

Here pnm is the spectral density of the resolvents, p„,„ — 

_ oo 

— -luiRnmii^) and Z = Tr J de e~^'^p{e) is the cluster 

— oo 

partition function. 

6 Results 

The DCA enables us to include short-ranged nonlocal dy- 
namical correlations neglected in the DMFA. The main 
goal of this section will be to show that this is indeed the 
case and to present some systematics on how these non- 
local correlations evolve and in what way their influence 
depends on system parameters like filling and band struc- 
ture. To this end we present results for the Hubbard model 
(|l|) on a square lattice in two dimensions. For nearest- 
neighbor hopping the dispersion is et = — 2i(cosfc2: -I- 
cosfcy), i.e. the bandwidth of the noninteracting system 
W = 8t. Calculations were performed for a 1x1 cluster 
{Nc = 1), which is equivalent to a DMFA calculation, and 
for a 2x2 cluster {Nc = 4) . A comparison of the results for 
both cluster sizes is used to study the eflFect of nonlocal 
correlations present in the Nc = 4, but neglected in the 
Nc = 1 calculation. 

The total number of cluster eigenstates scales with the 
cluster size Nc like 4''^'^ . The large number of eigenstates 
(256) for the 2x2 cluster results in an expensive numeri- 
cal calculation. The complexity of the problem can be re- 
duced by taking into account the symmetries of the cluster 
Hamiltonian (|l3|). Since our studies are restricted to the 
paramagnetic phase we can drop the spin index due to 



the SU(2) symmetry of the cluster Hamiltonian (tt3[). A 
further reduction in complexity can be achieved by using 
the point-group symmetry of the cluster. However, this 
depends strongly on the choice of the cluster momenta K 
within the first Brillouin zone. A priori, there is no re- 
striction in the choice of the cluster momenta K within 
the first Brillouin-zone, since in the derivation of the DCA 
algorithm no special assertion about the cluster K points 
was made. One e.g. could choose all K momenta to lie on 
the Fermi surface. However, to identify eigenstates which 
are degenerate due to the geometric symmetry, one has 
to classify the eigenstates according to the cluster mo- 
menta K. Since the cluster Hamiltonian ([l3| ) conserves 
the cluster momentum, its many-particle eigenstates can 
be classified according to their total momentum, which is 
just the sum of the momenta of the participating one par- 
ticle states. This approach restricts the freedom in choos- 
ing the cluster momenta K to exactly one possibility. The 
only set of cluster momenta K which form a group under 
addition is Kai = In, where I = 0,1 and a = a; or y. This 
set corresponds to periodic boundary conditions for the 
cluster in real space. With this choice of the cluster mo- 
menta we are able to classify the eigenstates according to 
their total particle number, total momentum, total spin 
and their z-component of the spin. The degeneracy in the 
cluster momentum points (0, tt) and (tt, 0) and the spin 
symmetry finally lead to an effective number of 123 non- 
degenerate eigenstates which have to be considered. Then 
effectively only resolvents with different energies occur in 
our calculations. 

The remaining numerical task of calculating the cou- 
pled equations (|2^) and (|2^) self-consistently becomes for- 
midable as the cluster size increases. Although the study 
of larger cluster sizes is in-principle possible [p3|, presently 
this restricts our calculations to a cluster size of Nc = 4. 
Also the evaluation of two-particle correlation functions 
is formally possible, but the associated numerical effort 
scales much worse with the cluster size Nc than calcula- 
tions on the one particle level. Hence our calculations are 
currently limited to one-particle Green functions. 

We will present results for local single-particle spectra 
as well as for the bandstructure. Since within the DCA 
we calculate the self-energy at the selected momenta K 
only, we need to perform a bilinear interpolation of the 
self energy between the cluster momenta K to calculate 
nonlocal spectra. We also show results for the self-energy 
at the Fermi-surface. The shape of the Fermi-surface is 
not a priori clear. In order to evaluate the Fermi surface 
we take the bilinear interpolated self-energy and calculate 
the occupation n(k) in momentum space along various 
directions in the Brillouin zone. The maximum value of 
|dn(k)/dk| along these directions then marks the Fermi- 
surface and we get the self-energy at the Fermi surface 
from the interpolated form. 

In the following we will concentrate at first on a generic 
set of values for the temperature T and Coulomb param- 
eter U, namely T = WjXh and U = W/2. These values 
for U and W assure that for the half filled case the sys- 
tem is metallic and far from the Mott-Hubbard transition, 
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which is expected to occur at values U Ki W. This choice 
allows us to directly compare the properties of the ex- 
pected metallic phases at and off half filling. The equally 
interesting question, of how the Mott-Hubbard transition 
at half filling will be affected by nonlocal correlations is 
left out for the time being and will be the subject of a 
forthcoming publication. 
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Fig. 6. Local density of states for the 15% doped 1x1 (N^ = 1) 
and 2x2 (A^c = 4) cluster at fixed temperature T = W/15 and 
interaction U — 0.5W . Inset: Imaginary part of the correspond- 
ing self energies in a narrow region around uj — 0. 



Fig. H shows the density of states for both the 1x1 and 
2x2 clusters for a doping oi S = 15%. Both spectra display 
qualitatively similar features, namely the typical Hubbard 
bands and an enhanced density of states at the Fermi level 
oj = 0. For both cluster sizes the imaginary part of the 
corresponding self energies shows a parabolic minimum 
at the Fermi level as expected for a Fermi liquid, where 
— Imi7(w)|^^g ex w^ -I- ttT^. In the case of the 2x2 clus- 
ter the self energy at the Fermi wave vector E(kF,t^) is 
obtained from the above mentioned bilinear interpolation. 
Here kp lies along the diagonal from (0, 0) to (tt, tt); how- 
ever, for this set of parameters the self energy on the Fermi 
surface is almost isotropic and does not change its qualita- 
tive behavior as a function of the wave vector on the Fermi 
surface. The influence of nonlocal short-ranged dynamical 
correlations is visible in the 2x2 cluster calculation only in 
a slightly enhanced scattering rate at the Fermi level and 
therefore in a slightly reduced density of states as com- 
pared to the 1x1 result. The additional structures on top 
of the Hubbard bands can be traced back to the complex 
multiplet structure of the cluster. However, the qualitative 
physics remains unchanged by increasing the cluster size. 
This observation is consistent with the fact that at such 
strong doping antiferromagnetic fluctuations have prac- 
tically died out and should thus show no significant in- 
fluence on the physics of the system. The appearance of 



the quasi-particle resonance at low enough temperatures 
is well known in the case of the 1x1 cluster (DMFA): There 
it was shown that the evolution of this quasi-particle res- 
onance with decreasing temperature is accompanied by a 
reduction of the effective local magnetic moment |p|,p2|. 
This interplay of both effects is a fingerprint of the Kondo 
effect occurring in the single impurity Anderson model, 
which underlies the DMFA. Our results thus suggest that 
for the lattice system the physics is quite similar and the 
quasi-particle resonance at the Fermi level reflects Kondo 
like physics. It is important to note that this means that 
the Kondo like behavior in the Hubbard model is, at least 
for moderately to strongly doped systems, a real feature 
of the model and not an artifact of the limit of large di- 
mensions. 

For weakly doped or half filled systems, short-ranged 
antiferromagnetic spin fluctuations will be present and 
strong even at temperatures well above a magnetic tran- 
sition. One thus expects that physics of the system will 
be strongly influenced and may even develop non-Fermi- 
liquid-like behavior. 




Fig. 7. Local density of states for the 1x1 {Nc = 1) and 2x2 
[Nc = 4) cluster at fixed temperature T = W/15, interaction 
U — 0.5W and half filling. Inset: Imaginary part of the corre- 
sponding self energies in a narrow region around lo — 0. 



Since these fluctuations will be strongest in the extreme 
case of half filling, we will consider this case next. Fig. |^ 
shows the results for both cluster sizes with the same pa- 
rameters as in Fig. ^ but at half filling. Whereas the 1x1 
cluster result displays the same features as in the doped 
case - enhanced density of states at the Fermi level accom- 
panied by a parabolic minimum in the imaginary part of 
the self energy - the spectrum of the 2x2 cluster calcula- 
tion is completely different from the doped case: Instead 
of forming a quasi-particle resonance as in the DMFA, the 
density of states develops a pseudogap at zero frequency 
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and the corresponding imaginary part of the self energy 
which is again almost isotropic displays a strongly en- 
hanced scattering rate at the Fermi energy. This surprising 
and interesting behavior has two possible explanations. 
The first and physically most appealing one is that short- 
ranged antiferromagnetic fluctuations do indeed drive the 
system from a Fermi liquid into a non Fermi liquid at tem- 
peratures high compared to the Neel temperature. Note 
that the underlying mechanism is very similar to the in- 
terpretation of the pseudogaps observed in the high-Tc 
compounds well above T^- The second interpretation is 
that the nonlocal corrections yield a reduction in the crit- 
ical value Uc at which the Mott-Hubbard metal-insulator 
transition occurs. 



U=0.5W 

T=W/15 
5=0 




Fermi-energy which is most pronounced at the X-point 
(7r/2, 7r/2). But in addition to these effects we now can re- 
solve additional incoherent background structures at the 
points r at E{k) k, 1 and M at E{k) k, -1. These addi- 
tional states are just shifted by the wave vector Q — (tt, tt) 
with respect to the main bandstructure. 



1.2 



k=(0,0) 

k=(7t,7t) 




Fig. 9. Spectral density at the points F and M for the same 
parameters as in Fig. [^ and y. Note the occurrence of shadow 
structures with small spectral weight in addition to the main 
structures. 



Fig. 8. Contour-plot of the spectral density for the 2x2 cluster 
calculation along the main symmetry directions as indicated in 
the inset for the same parameters as in Fig. M. The dark color 
marks regions with high spectral density. The open symbols 
indicate the maxima of the spectral function. The solid line 
marks the dispersion for the noninteracting system. 



A contour-plot of the spectral density ^(k, oj) obtained 
with the bilinear interpolation scheme discussed earlier 
along the main symmetry directions in the Brillouin zone 
is shown in Fig. pi The dark shading marks regions with 
high spectral density. The open symbols in Fig. @ repre- 
sent the positions of the most pronounced local maxima 
of A(k, w) and can be viewed as effective band structure 
of the 2D Hubbard model for the set of parameters under 
consideration. Compared to the bandstructure of the non- 
interacting system (illustrated by the solid line) the inter- 
actions have various effects on the spectrum. The band of 
the noninteracting system splits into two separated bands 
above and below the Fermi-energy. Note that the spectral 
density corresponding to the nearly dispersionless features 
of the two bands is very low and comparatively broad 
in the regions without states of the noninteracting sys- 
tem. We again notice the opening of the pseudogap at the 



In order to better resolve these structures we show the 
corresponding spectral density at the points F and M in 
Fig. p. In addition to the main structure at oj < for the 
point F {u) > Q for M) we notice satellites at w > for F 
(w < for M) with small spectral weight. These new states 
are absent in the non-interacting system as well as in 
the DMFA and result from the nonlocal antiferromagnetic 
correlations. Even in the paramagnetic phase the short- 
ranged antiferromagnetic spin fluctuations are sufficient to 
produce this indication of the ordered phase. Such a pre- 
cursor effect of the antiferromagnetic long range order can 
for example be seen in the cuprates ll^| in Fermi-surface 
measurements. The observation of these spin- fluctuation- 
induced shadow states accompanied by an opening of a 
pseudogap strongly supports the first suggested scenario 
of the antiferromagnetic spin fluctuations driving the sys- 
tem to a non Fermi liquid. 

To gain more insight in the nature of the pseudogap 
and elucidate the physics of the observed non Fermi liq- 
uid behavior we added a next nearest neighbor hopping 
t' to the hopping integrals tij in the Hamiltonian (nl). 
The dispersion then becomes ek — — 2t(cos k^ + cos ky) -\- 
At' cos kx cos ky and the Fermi surface is no longer nested 
because ek+q ^ — ek for the nesting vector Q — (tt, tt). By 
including t' we can thus frustrate the lattice and gradually 
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suppress antiferromagnetic spin fluctuations. On the other 
hand, since we keep the non-interacting bandwidth W and 
therefore the ratio U/W fixed when including the t', we 
do not expecte this change to affect the Mott-Hubbard 
transition very drastically. Therefore, if the pseudogap is 
a precursor of the Mott-Hubbard transition we do not ex- 
pect it to be influenced dramatically when we increase t' . 
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Fig. 10. Local density of states for the 2x2 (A'^c ~ 4) cluster at 
fixed temperature T — W/15, interaction U — 0.5W and half 
filling for various values of the next nearest neighbor hopping 
integral t' . 



interaction strength U — 1.3W and lower temperature 
T — W/30. For the hopping integrals we chose t — 0.25 
and t' — — 0.35f to imniitate the measured Fermisurface 
of underdoped cuprates in the non-interacting system [ |l2[ 
Q . Since we expect the spectra to display sharp features 
at this low temperature, we refrain from performing a bi- 
linear interpolation of the self-energy and show results for 
the coarse grained spectra only. Fig. ^ shows the coarse 
grained spectral functions A{K,llj) = — l/7r3mG'(K, w) 
at K = (0,0) and K = (tt, 0) for a doping 6 = 10% and 
temperature T = W^/30 in a narrow region around the 
Fermienergy. First we notice a strong anisotropy in the 



(O.ii) (n.K) 
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(7t,0) 



Fig. |l^ shows the behavior of the local density of states 
as the next nearest neighbor hopping t' and therefore the 
lattice frustration increases. First note that we lose par- 
ticle hole symmetry for t' ^ because the bare density 
of states is no longer particle hole symmetric. Obviously, 
as t' increases, the pronounced pseudogap for t' = gets 
gradually reduced. However, even for the maximum value 
I i'l = |i| we still see a small dip at zero frequency, which 
possibly can be attributed to the fact that even for a 
completely frustrated system extreme short-ranged fluc- 
tuations will be present, which are however strongly re- 
duced in magnitude as compared to the nested situation. 
We thus conclude that indeed the nonlocal, short-ranged 
antiferromagnetic spin correlations are responsible for the 
development of the pseudogap at the Fermi energy in the 
Hubbard model with t' = at half filling, which in fact 
must be viewed as precursor of an antiferromagnetically 
ordered state at much lower temperatures. One highly in- 
teresting question to be addressed in our future work will 
be of what precise nature this non Fermi liquid state is 
and how it might be related to several phenomenologi- 
cal scenarios proposed for the two-dimensional Hubbard 
model. 

As a further illustration of the ability of the DCA to 
include nonlocal correlations, we show results for a larger 



Fig. 11. Coarse grained spectral function of the cells around 
K = (0, 0) and (vr, 0) for Ac = 4 and a doping concentration 
S = 10%, temperature T = VK/30 and hopping integrals t = 
0.25 and t' = — 0.35i. The inset shows the Fermisurface of the 
corresponding non- interacting system (solid line), the nested 
k-points (dot-dashed line) and the boundaries of the coarse 
graining cells (dashed lines). 



coarse grained functions. For the states located in the cell 
around K = (0, 0) the spectrum is peaked at the Fermi- 
energy and only slightly broadened due to the finite tem- 
perature as characteristic for a Fermi liquid. The situation 
is completely different in the cell around the point (7r,0). 
All the spectral weight is transfered to broad features at 
higher energies, which represents the incoherent part of 
the spectrum. Therefore a pseudogap opens at the Fermi- 
energy. 

Our results are qualitatively in agreement with calcu- 
lations [p3 for the spin fluctuation model, in which the 
magnetic interaction between the quasiparticles is held 
responsible for the anomalous normal state properties of 
the cuprates. This method provides a direct explanation 
for the anisotropic behavior of the spectral density. Within 
this approach the electron-electron interaction is mediated 
by an imperially determined dynamical spin susceptibility. 
Since this susceptibility is strongly peaked at the antifer- 
romagnetic wave vector Q — (tt, tt), one has to distinguish 
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two different regions of the Fermi surface: Quasiparticles 
in regions of the Ferniisurface which can be connected by 
the wave vector Q are called hot quasiparticles because 
they feel the full effects of the spin fluctuation induced 
interaction because of nesting. This is illustrated in the 
inset of Fig. O, which displays the corresponding non- 
interacting Ferniisurface for the chosen parameters. One 
can notice that the hot quasiparticles are located in the 
cell around (tt, 0) and therefore represented by the coarse 
grained spectral function at K = (tt, 0). On the other hand 
the cold quasiparticles located along the diagonal couple 
only weakly to the spin excitations, since the Ferniisurface 
in this region is not nested. This part of the Femisurface 
falls in the cell around K = (0, 0) and therefore the spec- 
trum in this cell displays Fermi- liquid like behavior. The 
spectrum around (tt, 0) on the other hand gets strongly 
renormalized due to the strong coupling to the spin exci- 
tations. This phenomenological picture provides a direct 
explanation of our results qualitatively consistent with the 
calculations for the spin fluctuation model. 

Calculations for larger doping {6 = 20%) (not shown 
here) show that this effect of the anisotropic behavior of 
the spectrum and the opening of the pseudogap in the 
hot regions disappears. This observation can also be un- 
derstood within the picture of the spin fluctuation induced 
correlations, since the antifcrromagnctic spin fluctuations 
become strongly suppressed upon doping. 



opening of a pseudogap in the density of states and there- 
fore non Fermi liquid character. Calculations with a next 
nearest neighbor hopping t' show evidence that this pseu- 
dogap is due to antiferromagnetic spin correlations and 
therefore a single particle precursor of the antiferromag- 
netic phase transition. The pseudogap persists upon weak 
doping in qualitative agreement with the spin fluctuation 
model for the cuprates. 

Hence our calculations have shown that for the weakly 
doped system, nonlocal correlations play an important 
role on the single particle properties and change the char- 
acter of the system from a Fermi liquid to a non Fermi 
liquid. These non-local features are missing in the DMFA 
spectra, but appear in the DCA spectra as soon as the 
cluster size exceeds one. 
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7 Conclusion 

We motivated the recently introduced Dynamical cluster 
approximation (DCA) by its microscopic definition based 
on a choice for the Laue function. It partially restores the 
momentum conservation at the internal vertices which was 
relinquished in the Dynamical Mean Field Approximation 
(DMFA). The resulting theory maps the lattice problem 
onto a self-consistently embedded periodic cluster of size 
Nc- The DCA is a fully causal and systematic approxima- 
tion to the full lattice problem with corrections 0{1/Nc) 
in two dimensions. We develop a Non Crossing Approxi- 
mation (NCA) to solve the effective cluster problem which 
is a systematic ^-derivable approximation to the cluster 
problem with corrections 0{1/N^). 

We applied our DCA-NCA formalism to the Hubbard 
model on a square lattice and calculated the single-particle 
properties when Nc = 1 and Nc = 4. For a highly doped 
system, with S = 15% and a Hubbard U of half the bare 
bandwidth, nonlocal correlations present when Nc = 4 
turned out to be unimportant and both cluster sizes yield 
qualitatively similar results with an enhanced density of 
states and a minimum in the scattering rate at the Fermi 
energy. Thus, independent of the cluster size, the highly 
doped system showed Fermi liquid character. However, 
as the doping decreases, the non-local correlations be- 
come far more important, and the half- filled system dis- 
plays strikingly different results for the two cluster sizes. 
Whereas the Nc — 1 cluster result still displays Fermi liq- 
uid like behavior, the iV^ = 4 cluster calculation shows the 
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A On the analyticity of r{K, z) 

The following proof of the analyticity of r(K, z) is based 
on the derivation of an analytic expression for the cluster 
Green's function 



G(K,z) 



E' 



1 



N, 



N ^ z + ^ - EK+k' - ^c(K, z) 



The self-energy function Z'c(K, z) is assumed to be an- 
alytic in the upper and lower half of the complex plane 



with sign(3?7iZ'c(K, z)) = — sign3?n(z). Our proof em- 
ploys standard methods of projection technique |pq,B6[. 
To this end we abbreviate ^k = z+fjL—Sci)^, z), 3m(^K) > 
and introduce a N/Nc dimensional set of linearly in- 
dependent vectors |k') and a linear hermitian operator 
Hk satisfying Hxlk') = eK-fk'|k') and E|k')(k'| = 1. 

k' 

There obviously exists another vector |/o) with (k'|/o) — 
^jNc/N for all k'. With these conventions we may write 

1 



G(K,z) = (/o|- 



-|/o> 



The resolvent operator may be trivially rewritten as 



1 



1 



1 



TYk- 



1 



^K — Wk Ck Ck Ck — Hk 



and thus 
(/o' 



7 Ti—\h) = ^ + ^(/o|Hk7 :rr- 



l/o) 



where we made use of (/o|/o) = 1- We now define two 
projection operators P = |/o)(/o| and Q = 1 — P and 
insert \ = P + Q after TYk in the second term, leading to 

(/o L ^^ l/o) = ^ + ^(/o|Hk|/o)(/o L ^^ l/o) 
+^(/o|HkQ , ^^ |./o) . 

4k 4k — rlK 



Since furthermore 
1 



1 1 



Ck — Hk ^k — WkQ — 'HviP ^K — WkQ 



Ck — "^kQ 

and due to Ql/o) =0 

1 



^K — Hk 



' Ck — 'Hk.Q 



l/o) = 



we finally obtain 

(/oIt t;— l/o) = 7 ^' ■^(/o|'Hk|/o)(/q|7 TTT- 

4k — /Ik 4k 4k 4k — Axk 



l/o) 



+ ^(/o|^kQ7 nj ^ Hk|/o)(/o|7 TT 

4k 4k - 'XK^' 4k - /XK 

With Q^ — Q we may rewrite 

1 .... „ 1 



l/o) 



■HkQ- 



:H 



K 



^K — 'H-K.Q 

With the abbreviations 

r(K,z) = (/o|7^KQ 



n^Q 



Ck - QHkQ 



QHk 



1 



i^K — QTi-KQ 



O^kI/o) 
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and 



CK = (/oIWkI/o) = -j^ y ^ eK+k 



the final result is 



G(K,z) = 



U-^-K-r(K,z) 



(26) 



where we dropped the frequency argument for convenience. 
By defining 

^K+k = ^K+k ~ ^K , (29) 

such that -^c/^X^ic^K+k = we can make use of the 
exact relation 



Note however that the averaging procedure replaces the 
kinetic energy of the lattice e^ by a quantity coarse grained 
onto the cluster. 

We are now left with the proof that sign5r7i(/^(K, z)) — 
—sigDSm{z) is fulfilled. This however can easily be done 
by repeating the above step for the new vector |/i) = 
QHk|/o) = Hk|/o) - exl/o), (/o|/i) = 0, appearing in 
the definition of I^(K, z). Note that this is simply the first 
step in a Schmitt-Graham procedure to generate an or- 
thogonal set of vectors. With Hk — QT^-kQ it follows 



G(K + k) = g(K) + .g(K)tK+kG(K + k) 



(30) 



By inserting (|30| ) in (28) it is straightforward to show that 
r(K) is given by 



N 



E^K+kG(K + k) 



r(K) 






(31) 



kG(K 



r(K,z) 



bl 



•?K - fli - r{K, z) 



(27) 



where 



ai 



>0 



(/il/i) 
(/ol/o) 

(/i|Hk|/i) 



By performing a Taylor series expansion of t^,^ around 
the cluster points K it can be seen that tj^,^ is of or- 
der 0(Ak), where Ak = I-kJNJ . Therefore we see from 
Eq. H^that r{K) is of order O^Akf) and for two dimen- 
sions we finally get the result 



r(K) 



0(\IN,) 



(32) 



f(K,z) 



(/il/i 



-(/iIHkQ 



1 



Ck - QHkO 



-Q^k|/i 



where Q now projects onto the subspace orthogonal to 
l/o) and l/i). It is clear from the above result that this 
procedure can be repeated, leading to a sequence of mu- 
tually orthogonal vectors |/„) = 7iK|/«-i) -an-i|/n-i) - 
6^_2|/n-2) and acontinued fraction representation of G(K, 
with coefficients 



^) 



^Li = 



\Jn\jn) 



(/n-l|/n-l) 



>0 



_ (/nl^Kl/n) 
""" ifnlfn) 

for n > 1. It is important to emphasize that the re- 
sulting coefficients &^„i are non- negative by construction. 
This however immediately leads to the desired relation 
sign5TO(r'(K, z)) — — sign577i(z) and hence to the causal- 
ity of r(K,z). 



B Proof of r(K) ~ 0{1/N^) 

The following proof is based on the definitions (|l^JTl|) of 
the coarse grained Green function G. With g^^{K.) = uj — 
ck — ^(K) we can rewrite G in the form 



.g-i(K)-r(K) 



(28) 



